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Abstract 



It is shown, contrary to common belief, that the Fritzsch ansatz 
for the quark mass matrices admits a heavy top quark. With the 
ansatz prescribed at the supersymmetric grand unified (GUT) scale, 
rS I one finds that the top quark may be as heavy as 145 GeV, provided 

j^ ' that tan/? (the ratio of the vacuum expectation values of the two higgs 

doublets) S> 1. Within a non-supersymmetric GUT framework with 
two (one) light higgs doublets, the corresponding approximate upper 
bound on the top mass is 120 (90) GeV. Our results are based on a 
general one-loop renormalization group analysis of the quark masses 
and mixing angles and are readily applied to alternative mass matrix 
ansatze. 
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The standard SU{3) x SU{2) x f/(l) gauge model predicts the existence of 
a sixth (top) quark whose mass is known experimentally to be heavier than 
91 GeV.^ Analysis based on electroweak radiative corrections may favor a 
moderately heavy top quark, say below 160 GeV, with a central value possibly 
around the 120-140 GeV mark.^ Other independent sources of information, 
such as recovering the 6 quark mass in the range 4.25 ± 0.1 GeV within 
a supersymmetric grand unification framework (SUSY GUTs) with unified 
third generation Yukawa couplings, also lead to a top mass in the above 
range. ^ 

A simple approach for incorporating the observed hierarchy of quark 
masses and their mixings was suggested a long time ago by Fritzsch.^ It 
prescribes a form for the mass matrices which has a certain amount of pre- 
dictive power. Based on considerations of the Vcb element of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, it has been argued that if the top quark 
mass turns out to be significantly heavier than 90 GeV, then the simplest 
version of the Fritzsch ansatz is excluded. ^'^ 

The main purpose of this note is to point out that the Fritzsch ansatz 
actually permits the top quark to be much heavier, provided it is considered 
within the framework of grand unification involving more than one 'light' 
higgs doublet. Perhaps the most prominent example of this is provided by 
supersymmetric grand unification. An important new parameter is tan/?, 
the well known ratio of the vacuum expectation values of the two higgs dou- 
blets that supersymmetry mandates. It turns out that with tan/3 sufficiently 
greater than unity, the grand unified supersymmetric version of the Fritzsch 
ansatz permits the top quark to be as heavy as 145 GeV. A similar result 



holds within a non-supersymmetric GUT framework provided that there are 
two 'hght' higgs doublets. In this case the top quark mass can be as high 
as 120 GeV. For completeness, we also present the results for a single higgs 
doublet and confirm that the top mass cannot exceed by much 90 GeV. 

Our analysis is based on the one-loop renormalization group evolution 
of the quark masses and mixing angles. The results are quite general and 
can be used to test alternative mass matrix ansatze as well. As an additional 
application, we derive lower limits on the top mass in the standard, two-higgs 
and supersymmetric models assuming that the relation \V^,j\ = Jm^/m^ 
holds at the unification scale. 

The Fritzsh ansatz for the quark and charged lepton mass matrices takes 
the following form at the GUT (or some other appropriate superheavy) scale: 

/ au,d,e \ 
Mu4,i = Pu,d,e 0'u,d,e bu,d,e Qu,d,e (1) 

\ bu,d/ Cu,d/ J 

where Pu,d,e and Qu,d,e denote diagonal phase matrices and a,b,c are real 
(positive) quantities. Among the various phases contained in P and Q, only 
two are relevant for quark mixing, we denote them by ip and 0. The ansatz 
predicts the CKM matrix elements in terms of the quark mass ratios and 
these two phase parameters. In order to incorporate the fermion mass hier- 
archy, one finds that c ^ b ^ a. From the mass eigenvalues in the quark 
sector in particular, 

Cu — rn^, bu — (m°m°)2, a„ ~ (m°m°)2 
Cd — m^i bd — {m^m'l) "2 , ad — im^m^)^ 



The superscripts are to emphasize that the ansatz (1) is prescribed at some 
superheavy scale. [As a consequence of grand unification the parameters 
o-d-, bd, Cd of the down sector typically are related to a^, b^, q of the lepton 
sector.] 

Our main concern here is the following asymptotic relation which is a 
consequence of eq. (1): 
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(3) 



The implications of (3) at low energies are evaluated by studying the evolu- 
tions of the gauge and Yukawa couplings. In the absence of Yukawa couplings 
that are comparable to or larger than the gauge couplings, relation (3) would 
be essentially unrenormalized, and would require that the top mass not ex- 
ceed 90 GeV. However, since the top quark mass is known to be greater than 
91 GeV, and indeed may be significantly larger than this, its Yukawa cou- 
pling cannot be ignored. Moreover, models in which the up and down type 
fermions obtain masses from couplings to distinct higgs doublets allow for 
large intrinsic Yukawa couplings hi, and h^-. We will exploit, in particular, the 
fact that if tan (3 (the ratio of the two doublet vev's which provide masses to 
up type and down type fermions) is sufficiently large, then /ib(/ir) can even 
exceed /ij! As a consequence, in contrast to the single light higgs doublet 
case, relation (3) it turns out can lead to a low energy prediction for | Vcb \ 
which is compatible with the existence of a heavy top quark. 

The Yukawa sector relevant for our discussion has the generic form 



£y = qLHu(l)uUR + qLHd(j)ddR + iLHi(f)deR + h.c. 



(4) 



where H^, Hd, Hi denote the 3x3 Yukawa couphng matrices for the up 
quarks, down quarks and charged leptons. Recall that in the standard model, 
0u = (pd- The one loop evolution equations for the Yukawa matrices take the 
form (t = in{fi/Mz)y.^ 
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+ -{bHdHl + cH.,Hl)-GD 



Hd 
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(5) 



For the three low energy models under discussion the coefficients a, 6, c are 
given by 



(a, 6, c) 



(0'2,i; 



MSSM (Minimal SUSY) 

Two higgs doublets (non — SUSY) v"/ 

SM (Standard Model) 



= (0,1, i) 

= (1,1,-1) 
Also, the quantities G^/, Gd and Ge for the SUSY (non-SUSY) case are 

respectively: 



Gu -- 


- f^? + 3^i + f^3^ 


Gd -- 


= l9l + ^9l+'i9l 


Ge -- 


= ^9l + ^9i ; 



(MSSM) (7) 



Gu = M + l92 + 89i 

Gu = M + 19I + 89I (Two higgs ; SM) (8) 

^E — —9l + 492 ■ 

The gauge couplings Qi (above) obey the standard one loop renormaliza- 



tion group equations: 



^^''^ = h.9t. ^ = 1,2,3 (9) 



where 



(h^MM) = (11,1,-3) MSSM 

= (7,3,-7) Two higgs (10) 

= (f,^,-7) SM 
Note that we are not using the S'f/(5) normalization for the hypercharge 
coupling Qi. From eq. (5), one can compute the evolution equations for the 
eigenvalues of the Yukawa coupling matrices: ^'^ 



i« 2% 
lovr — — 

dt 



3 Y. 9] + ^a J2 9l + a ^ 9b ' Gu 

j=u,c,t l3=d,s,b b=e,n,T 



IGtt 



2^9a 



dt 



l3=d,s,b 



3a E ^J + 3 E 9l+ E 9b-G 

j=u,c,t f3=d,s,b b=e,fi,T 
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■i^bgl + T^c E ^1 I V, 



ja 



j=u,c,t 



IQtt' 



,dga 
dt 



9a 



3a E ^1 + 3 E ^?+ E 9'b-G, 

j=u,c,t P=d,s,b b=e,fi,T 



3, 2 

+ -,b9l^ 



(11) 



where i = {u,c,t), a = {d, s,b), a = (e, /i, r). 

We will also need the evolution equations for the elements of the CKM 
matrix:'^'^ 



167r^- I V, 
dt 
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■ (12) 



The above expressions simplify considerably if we exploit the hierarchy in 
the Yukawa couplings {gb ^ Qs ^ gd, etc) and in the CKM matrix elements. 
If only the leading terms are kept, one obtains the following approximate 
expressions for the evolution of the various mass ratios and the mixing angles: 
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Several comments are in order: 

1. The quantities | \4b |, | Kb I, I ^td I and | Vts \ have identical evolutions. 

2. The evolutions of | Vus \ and | Vcd \ involve second family Yukawa 
couplings and consequently are negligible. 

3. iTLd/mi, and ms/mf, have identical evolutions. Similarly for rriu/mt and 
rric/mt. 

4. The evolutions of md/ms and mu/rric can be ignored. 

5. The CP violating parameter J is not independent, its evolution can be 
obtained from the running of the square of |V"cb|, for example. 

Using eq. (13) the asymptotic expression (3) can be recast in the form 



Vcb I V, 



cb 



VsbJ e^VtW — 

V rrih V TTit 



(14) 



where \Vcb\ stands for the CKM angle at the weak scale, and the r^'s are 
the respective running factors, rjcb = lKSl/|Kfe|, Vsb = {^)/{t^), etc. To 
compute these rj factors, one has to solve numerically the evolution equations 



in eq. (13). The running of gt,gb,gT are obtained from eq. (11), neglecting 
the hght two famihes. As input parameters we take 



(15) 



For the superheavy scale Mx above which the asymptotic relations are valid 
we take 



ai{Mz) -- 


= 0.01013 


a2{Mz) = 


= 0.03322 


ocz{Mz) = 


= 0.115 . 



M 



X 



(16) 



10^6 GeV (SUSY) 

= 10^4 GeV (non - SUSY) . 
For the running quark masses we use the values of ref. (9). In particular, 
m,(l GeV) = 175 ± 55 MeV, mjym^ = 1.27 ± 0.05 GeV and mb{mb) = 
4.25 ±0.1 GeV. 

To get a feeling for the effect running has on the prediction for \Vcb\, let 
us first define the functions 

1 fln{Mx/Mz) 



fth = exp 



9t,bi^)dT 



(17) 



IGtt^j 

Note that ft^b ^ 1- Let us first consider the supersymmetric case. For small 
tan/3, Qb can be neglected since gb ^ gt- The renormalized prediction for 
\Vcb\ in terms of the mass ratios at the weak scale is then 



w 



cb\ 



nib V ^t 



(18) 



Since ft > ^, "v^e see that the disagreement of \Vcb\ with observations is more 
prominent for small tan/3. This feature persists as tan/3 is increased, until it 
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approaches mt/mi,. For tan/? = mt/mi), gt ~ Qh at all scales since both cou- 
plings track the same evolution equations (except for small corrections from 
the hypercharge and the r Yukawa couplings). To a good approximation, we 
see that the renormalized |V"cb| obeys 



IK 



cb\ 




m, _ ^i^ im, 




(19) 



That is, the relation (3) is essentially unrenormalized and the top quark 
cannot be heavy. This prompts us to consider values of tan/3 even larger 
than mt/mi), in which case better agreement is possible. Indeed, if we take 
the extreme limit gt ^ Qt, then in terms of fb defined in eq. (17), we see that 



Wr 



cb\ 



y rrib y mt 



(20) 



Note that the first term in eq. (20) has a reduction factor whereas the second 
term is enhanced, precisely what one needs to bring the prediction for \Vcb\ 
down and allow for a heavy top quark. These features are indeed borne out 
by the actual numerical computation to which we now turn. 

In fig. la, we plot the behavior of the ratios gt{M x) / gt{,M z) , gb{M x) / gb{,M z) 
and gr{Mx)/gT{Mz) as functions of rrit for tan/? = 3. Throughout this pa- 
per, we shall mean by rrit the running mass mt{mt). The pole mass is related 
to the running mass via 



m 



pole 



mt{mt) 



4 aq 

1 + 7T- 
3 TT . 



(21) 



Fig. lb shows the dependence of the functions f{Mx)/f{Mz) for / = 
iKfel) liTT^s/rribl and \mc/mt\. The behavior of \Vub\, \Vtd\ and \Vts\ are iden- 
tical to that of \Vcb\- Similarly, \mii/'mb\ runs as \ms/mb\ and \mu/mt\ as 
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\mc/mt\. The running factor for the CP parameter J is obtained from the 
square of \Vcb\- Since tan/? is small, the evolution of \ms/'mi)\ coincides with 
that of \Vcb\ which can be seen from eq. (13). In figures 2a and 2b, we plot 
the same quantities for the case of tan/? = mt/mb. Note that the running 
of Qt and Qb are almost identical. Finally, in figures 3a and 3b we plot these 
functions for tan/9 = 60. This value of tan/? corresponds to the infra-red 
fixed point solution for gb and Qr- 

From figures 1-3, one can compute the running factors rj^s entering in 
(14). It is clear from the discussions above as well as from figures 1-3 that 
unless tan/9 is larger than mt/mb, the situation for \Vcb\ is worse than the 
case with no running. However, for sufficiently large tan/9, |Vc6| is in the 
experimentally allowed range^*^ of 0.034 < \Vcb\ < 0.054. We give in Table 1 
the respective running factors for |Vc6|, {mg/mb) and {nic/mt) for tan/9 = 60 
and also list the lowest allowed value of | Vcb|. To arrive at the latter, we choose 
m,(l GeV) = 120 MeV, mb{mb) = 4.35 GeV, m^m^) = 1.32 GeV and set 
the phase = 0. We use two-loop QCD renormalization group equations 
to extrapolate the light quark masses from low energies to Mz- There is an 
upper limit of about 145 GeV on the top mass in this case coming from the 
requirement that Qb and Qr should remain perturbatively small in the entire 
range from Mz to Mx- One sees that for all values of m,t up to 145 GeV, 
\Vcb\ is in the experimentally allowed range. 
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rrit 


Vcb 


Vsb 


Vet 


1 '^cb |min 


80 


0.790 


0.510 


0.763 


0.023 


90 


0.781 


0.499 


0.746 


0.029 


100 


0.770 


0.485 


0.727 


0.034 


110 


0.757 


0.467 


0.703 


0.038 


120 


0.740 


0.445 


0.675 


0.042 


130 


0.718 


0.415 


0.639 


0.045 


140 


0.684 


0.371 


0.591 


0.046 


145 


0.660 


0.338 


0.559 


0.045 



Table 1. 

The running factors for |Vc6|, {ms/mb) and {mc/mt) in the 

super symmetric model with tan/? = 60 as functions of mt- 

In the last column is listed the lowest allowed value of |V"cfe|. 



The precise values of mt, \Vcb\ etc. will depend somewhat on the input 
value of as{Mz). For larger values of 0:3, the top quark can be heavier by a 
few GeV. For e.g., with as{Mz) = 0.12, the top mass is raised by about 5 
GeV. Most of the dependence arises from changes in extrapolating the light 
quark masses from 1 GeV to Mz- 

In figures 4,5 and 6, we plot the same quantities for the two higgs model 
for tan/? = {Sjint/mi,, 70). The behavior of the various functions is identical 
to the case of SUSY, except that due to the smaller beta function coefficients, 
the variation is not as pronounced. In Table 2, we list the running factors 
corresponding to tan/5 = 70 (the infra-red fixed point for gb and Qr) and 
conclude that mt as large as 120 GeV is allowed. 
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rrit 


Vcb 


Vsb 


Vet 


1 '^cb |min 


80 


0.887 


0.708 


0.875 


0.034 


90 


0.885 


0.705 


0.869 


0.040 


100 


0.882 


0.701 


0.862 


0.046 


110 


0.878 


0.697 


0.854 


0.050 


120 


0.874 


0.692 


0.844 


0.054 


130 


0.870 


0.686 


0.834 


0.058 



Table 2. 

Running factors and lowest allowed \Vcb\ in the 
two higgs model with tan/3 = 70. 

For completeness, we also display in figures 7a and 7b the variation of 
the relevant quantities in the standard model. Table 3 shows the r] factors 
from which it is clear that rrit cannot much exceed 90 GeV. This feature can 
also be understood qualitatively in terms of the function ft of eq. (17). The 
renormalized relation for \Vcb\ in the standard model case is 



IK 



cb\ 



f, 



rUb V iT^t 



(22) 



Although there is an overall suppression factor, the second term also becomes 
small and so the cancellation between the two terms becomes inefficient. 
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rrit 


Vcb 


Vsb 


Vet 


1 '^cb |min 


80 


1.020 


1.020 


0.981 


0.049 


90 


1.026 


1.026 


0.975 


0.054 


100 


1.033 


1.032 


0.969 


0.059 


110 


1.040 


1.040 


0.961 


0.064 


120 


1.049 


1.049 


0.953 


0.067 



Table 3. 

Running factors and |Vc6|min for 
the case of standard model. 



Turning now to some other predictions of the Fritzsch ansatz, the CKM 
elements \Vus\, \Vub\ and \Vtd\ have the asymptotic forms 



IV 



I ubl 



td 




V'\ = -^J^ + e'^, ^U -^~e'^, ^ 



(23) 



The first relation involves only the first two family masses and is therefore 
essentially unrenormalized. The phase ip should be near 7r/2 for agreement 
with the Cabibbo angle. From the second relation in (23), one can easily 
write down the renormalized value of \Vub\- It turns out that for tan/? = 60 
and rrit = 130 GeV in the SUSY model, the magnitude of the first term is 
only about 15% of the second term. Given that the phase ip ^ 7r/2 to a good 
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approximation, the first term can be neglected, which leads to the weak scale 

relation 

^~ /^~0.06 , (24) 

\Vcb\ V "^c 

in good agreement with the recent charmless B decay data. The same con- 
clusion also follows for the two higgs doublet model. As for \Vtd\, again one 
finds that the first term is negligible compared to the second term, so that 
at the weak scale we have 

iVtdl [md , s 

■wr^-\ — -0-2 • 25) 

\Vcb\ V "^« 

Finally, the renormalized value of the CP violating parameter J at the weak 
scale is given by (for ip = 7r/2, = 0) 



J ^ UhL Hhiiv^^l^ ^ S y, 10-^ ^ (26) 

consistent with observations. 

In the lepton sector, the Fritzsch ansatz can lead to two successful mass 
predictions provided that ad = ai, Cd = Q in eq. (1). Such relations arise 
naturally in GUT's if the elements a, c arise from a higgs 5 of SU{5) or a 
higgs 10 of 5*0(10). Note that bd and bi should be independent. The two 
asymptotic mass relations are 

m^ — ?T2g + m^ = 772° — m^ + m'^ ; m'^rri'l'm^ = m^m^rn^ . (27) 

The first relation would lead to a prediction of mii{mi,) = 4.2 GeV in the 
SUSY model with rrit = 130 GeV and tan/9 = 60, which is in good agree- 
ment with the value derived in ref. (9). The second relation would lead to 
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171^(1 GeV) ~ 7 MeV (if itls — 140 MeV), also in good agreement with 
observations. 

Since our analysis has been quite general, it is readily applicable to other 
mass ansatze as well. One particularly interesting case which has recently 
attracted a fair amount of attention is the asymptotic relation^^'^^'^^ 

(28) 




The renormalized relation in this case can be written down as 



ricb\Vc 



cb\ 



h hnc 
Vct\ — 



(29) 



In Table 4, we list rjcb, Tjct and the corresponding renormalized values of \Vcb\ 
for the case of the standard model. It is clear that there is a lower limit of 
about 170 GeV on the top quark mass in this case. 



rrit 


Vcb 


Vet 


\Vcb\ 


100 


1.033 


0.969 


0.081 


130 


1.060 


0.944 


0.068 


160 


1.104 


0.906 


0.058 


170 


1.124 


0.890 


0.054 


180 


1.149 


0.870 


0.051 


190 


1.180 


0.847 


0.048 


200 


1.221 


0.819 


0.044 



Table 4. 

Test of the asymptotic relation |\/°| = Jm^^/m^ in the 
standard model. mJirir) = 1.32 GeV is used. 
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How does this asymptotic relation fare in the supersymmetric and two 
higgs models? In the SUSY model, the value of \Vcb\ at weak scale can be 
written down as 

\Vc>,\ = f;h'\[^ . (30) 



f-2 f2 r^c 
Jt Jb \ • 

V ^t 



Since both ft and ft, are greater than unity, it is clear that smaller values 
of tan/3 will give better agreement. A similar conclusion holds for the two 
higgs model as well. In Tables 5 and 6, we list the r^'s and values of \Vcb\ as 
functions of rrit in the SUSY model and in the two higgs model for tan/5 = 3. 
In the SUSY model, the top quark mass should be close to the infra-red fixed 
point value of about 185 GeV. In the two higgs model, we see that there is 
no acceptable solution with rrit < 200 GeV and tan/5 = 3. Larger values of 
tan/5 are not displayed, since the agreement is worse. 



mt 


Vcb 


Vet 


\Vcb\ 


100 


0.972 


0.919 


0.084 


130 


0.946 


0.847 


0.072 


160 


0.895 


0.717 


0.063 


170 


0.862 


0.641 


0.060 


180 


0.800 


0.511 


0.057 


185 


0.708 


0.356 


0.052 



Table 5. 

Test of the asymptotic relation |Vj,^| = Jmf^/m^ 
in the SUSY model with tan/5 = 3. 
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rrit 


Vch 


Vet 


\Vct\ 


100 


0.988 


0.965 


0.084 


130 


0.978 


0.936 


0.074 


160 


0.963 


0.892 


0.066 


170 


0.955 


0.872 


0.064 


180 


0.947 


0.848 


0.061 


190 


0.935 


0.818 


0.059 


200 


0.921 


0.780 


0.057 



Table 6. 

Test of the asymptotic relation \V^ 



m^/m^ 



in the two higgs model with tan/3 = 3. 



To summarize, the Fritzsch ansatz for the quark mass matrices permits 
the top quark to be in the mass range suggested by the recent analysis of pre- 
cision electroweak data as well as by other independent estimates. However, 
it requires that the parameter tan /3 be considerably greater than unity, even 
exceeding mt/mh. It would be interesting to see whether this can be recon- 
ciled with the scenario of radiative electroweak breaking in the SUSY model 
which usually requires that tan/? be < mt/mi,. The ansatz also predicts 
values for the CKM matrix elements \Vus\, \Vub\, \Vtd\ and the CP violating 
parameter J that are in agreement with observations. 
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Figure Captions 

Fig. la: Plot oigt{Mx)/gt{Mz) (solid), gb{Mx)/ gb{Mz) (dashed) and 
gr{Mx)/ griMz) (dot-dash) as functions oi mt{mt) in the 
super symmetric model with tan/3 = 3. 

Fig. lb: The running factors f{Mx)/f{Mz) for / = \Vcb\ (solid), \ms/mb\ 
(dashed) and \mc/mt\ (dot-dash) versus rrit. The running factors 
for \Vub\, \Vtd\ and \Vts\ are identical to that of \Vcb\- Similarly, 
liTT'd/iTT'bl runs as Ims/mhl and \mu/mt\ as \mc/mt\. 

Fig. 2a: Same as in fig. la, except that tan/9 = mt/mh. 

Fig. 2b: Same as in fig. lb, except that tan/? = mt/mb. 

Fig. 3a: Same as in fig. la, except that tan/? = 60. 

Fig. 3b: Same as in fig. lb, except that tan/? = 60. 

Fig. 4a: The running factors for the two higgs (non-SUSY) model with 
tan/? = 3. Notation same as in fig. la. 

Fig. 4b: Running factors for the two higgs (non-SUSY) model with tan/? = 3, 
notation same as in fig. lb. 

Fig. 5a: Same as in fig. 4a, except that tan/? = rrit/mh. 

Fig. 5b: Same as in fig. 4b, except that tan/? = mt/mb. 

Fig. 6a: Same as in fig. 4a, except that tan/? = 70. 

Fig. 6b: Same as in fig. 4b, except that tan/? = 70. 

Fig. 7a: The running factors for the minimal standard model in the same 
notation as fig. la. 

Fig. 7b: Running of the mixing angles and mass ratios in the standard 
model, notation same as in fig. lb. 
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